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Abstract 

A simple two-level model is developed and used to test the proper- 
ties of effective interactions for performing nuclear structure calcula- 
tions in truncated model spaces. It is shown that the effective many- 
body interactions sensitively depend on the choice of the single-particle 
basis and they appear to be minimized when a self-consistent Hartree- 
Fock basis is used. 



1 Introduction 



Assuming that only two-body interactions act among nucleons in nuclei, we 
can write the nuclear Hamiltonian as 

i=i i<j 

This Hamiltonian can be divided into two terms, Hq and Hi, as 

ff = E(f^ + ^) + (E%--I>) =H + Hj, (2) 

i=l \ / \i<j i=l ) 

where Hq is the one-body Hamiltonian which defines a single-particle basis 
and Hi is the two-body "residual" interaction. If the Schrodinger equation 

H%(l,2,...,A)=Ei^i(l,2,...,A) (3) 

could be solved in the infinite Hilbert space for a many-body system, the 
results would not depend on the choice of the one-body potential u, or, 
equivalently, the single-particle basis. In practice, in order to solve Eq.(^), 
one must truncate the infinite Hilbert space to a finite model space and 
introduce an effective interaction (V e a) to be used in the truncated model 
space. For an yl-nucleon system, the effective interaction V e s will, in princi- 
ple, have two- [V^], three- [V^], and A-body [V^^] parts. If the exact 
effective interaction, containing all (two- to j4-body) components, can be ob- 
tained, the results will again be independent of the choice of the one-body 
potential u. However, three- and more-body effective interactions are diffi- 
cult to calculate and are often ignored in practical shell-model calculations 
with the hope that they are small. One is then left with only the two-body 

(2) 

effective interaction for performing nuclear structure calculations in a 
severely truncated model space. 

With recent developments on the effective interaction theory, i.e., the 
use of the no-core model space |[|, Q and the Lee-Suzuki approach H, |], ||] 
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to the folded-diagram series, we can now calculate the effective two-body 
interaction accurately. The question is how to choose u to give the best 

(2) 

approximation to the exact results when retaining only the V^a part of the 
effective Hamiltonian. 

In this work, we will use an exactly soluble two-level model (so that 
exact eigenenergies can be obtained) to show that the contribution to the 
ground-state energy from the often-neglected many-body effective interac- 
tions depend quite sensitively on the choice of the one-body potential u. We 
show that the Hartree-Fock (HF) self-consistent one-body field, appears to 
minimize the three- and more-body effects and consequently, the two-body 
effective interaction alone becomes a good approximation. 

This paper is organized as follows. In the next section, we give a brief re- 
view of the effective interaction theories. In section 3, we introduce a simple 
two-level model and calculate the energy-independent two-body effective in- 
teraction for an arbitrarily chosen, non self-consistent, single-particle basis. 
In section 4, we repeat the calculation performed in section 3 for a self- 
consistent Hartree-Fock basis and compare the results with those obtained 
in section 3. Finally in section 5, we give our conclusions. 



2 Energy-Independent Effective Interactions 

When the full Hilbert space is divided into P and Q spaces with P projecting 
out the model space from the full Hilbert space and Q projecting out the 
excluded space, the P-space effective interaction for an eigenstate with an 
eigenenergy E can be written || 

V cS (E) = PHjP + PH I Q E _ 1 QHQ QH I P , (4) 

Note that E, the energy of the many-body system which we wish to calcu- 
late, appears in the right-hand-side of the above equation. For this reason, 
the above effective interaction is also referred to as the energy-dependent 
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effective interaction. 

About ten years ago, Lee and Suzuki |J proposed a method for calcu- 
lating an energy-independent effective interaction. This method is based on 
a "generalized" G matrix, which is defined as 

G(co) = PHjP + PHjQ L^QHjP . (5) 

io — QHQ 

This generalized G matrix, which depends on a starting energy u>, is often 
referred to as the Q-box |7]]. Obviously, the generalized G matrix, or the 
Q-box, becomes the energy-dependent effective interaction V e g(E) when to 
is chosen to be equal to E. 

It is easy to see that for a two-particle system (A=2), the generalized G 
matrix becomes the Brueckner two-body G matrix Exact methods || [[(]] 
exist for calculating the Brueckner G matrix as a function of the starting 
energy u. However, the generalized G matrix is an A-body operator and is 
generally difficult to evaluate when A > 2. If it ever becomes possible to 
calculate G(oj), the energy-independent effective interaction V e g can then be 
obtained through an iterative procedure proposed in Ref.Q. Applying the 
Lee-Suzuki iterative procedure is equivalent to summing over all the folded 
diagrams which are not included in the generalized G matrix calculation. 
Further details on the Lee-Suzuki iterative methods can be found in Refs.||, 

II- 

In the case of a one-dimensional model space (one can always choose the 
model space to be one-dimensional no matter how many particles are in- 
volved), the effective interaction V e s for the model space becomes a number 
(i.e., a one- by-one matrix). It is then apparent by comparing Eq.(|j) and 
Eq. (|5|) that when the starting energy u> happens to be one of the eigenener- 
gies of the system, E, one has 

Veff = V eS (E) = G(uj = E) , (6) 

i.e., there is no need to utilize the iterative procedure and, therefore, the 
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contribution from the folded diagrams vanishes. In this case, no iterations 
are needed to obtain the effective interaction V e s from G(oS). 

In other words, in the case of a one-dimensional model space, the Lee- 
Suzuki approach to V e s simplifies to the following equation for E: 

E = PH P + G(E) , (7) 

which may be solved graphically. This has been shown explicitly in Ref.||]. 
It is also shown in Ref.||] that Eq.(^) is equivalent to the secular equation 
for the eigenenergies: 

det(il — EI) = . (8) 

Therefore, the solutions to Eq.(|?]) correspond to the exact eigenenergies 
of the A-body system. Where it has been tested the method appears to 
work well. Of course, it can only produce eigenvalues for states having 
non-vanishing overlap with the model-space state. 

Note that when the model space is .D-dimensional (D > 1), the energy- 
dependent effective interaction V e g(E) is different from the the LS energy- 
independent effective interaction V^ff . Both the model space Hamiltonians 
[Hq + F e ff (i?)] and [Hq + V e g] have D eigenvalues. But for the former Hamil- 
tonian, only one [E, i.e., the one at which V e a(E) is evaluated] out of D 
eigenvalues corresponds to a true eigenenergy of the system while, for the 
latter Hamiltonian, all its D eigenvalues are true eigenenergies of the system. 

Because it is not known how to compute G(co) for a system of A nucle- 
ons, the usual procedure is to start with the nuclear (or Brueckner) two-body 
reaction matrix G^ S and try to calculate the two-body effective inter- 

(2) 

action IQj . Assuming that the effective three- and higher-body forces are 
small, one can then use such an effective two-body interaction as input to 
shell- model programs, such as the OXBASH code [11] to perform model-space 



diagonalizations for the entire system. We briefly outline how one obtains 
the effective two-body interaction V^ s from the starting-energy-dependent 
Brueckner two-body G matrix. 
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The Brueckner G matrix represents the infinite summation (ladder sum) 
of two-particle scatterings. This is defined, in analogy to Eq.(|5|) for G(uj), 
as 

G^ (lj 2 ) = P 2 H I P 2 + P 2 H I Q 2 * — Q 2 HjP 2 , (9) 

lo 2 - Q 2 HQ 2 

where P 2 and Q 2 are now two-particle projection operators, which determine 
the allowed and forbidden intermediate states into which the two particles 
can scatter. Note that G^ is defined to depend on a starting energy u> 2 . It 
is clear that this starting energy refers to a two-particle system and when 
one approximates V e s by omitting many-body interactions, to 2 is different 
from the starting energy uj in G(co). 

(2) 

The ^-independent two-body effective interaction V^ s can be obtained 
from G^{oj 2 ) by summing over all the folded diagrams with two valence 
lines. This can now be accomplished without much difficulty by applying 
the Lee-Suzuki iterative procedure || || . 

3 v}g for an Arbitrary Basis 

A simple two- level model, consisiting of two single-particle states |1) and 
1 2), is used for our investigation. Each level can hold up to four nucleons: 
spin-up proton (p |), spin-down proton (p |) spin-up neutron (n f) and 
spin-down neutron (n j). The Hamiltonian for this model is determined by 
Hq, represented by the single-particle energies (spe) of the two levels (ei 
and e 2 ), and Hi, represented by 14 antisymmetrized, normalized two-body 
matrix elements (tbme) [(ab\Hi\cd) j^t]- Our choice for the SPE and the 
tbme is given in Table 1. 

As our first calculation, we diagonalize, using the oxbash shell-model 
code [JO]], the Hamiltonian in the full two-level space for A=2, 3 and 4. The 
results obtained are exact and are given in Table 2 in the column under the 
heading "Exact". 
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Table 1: The one-body part Hq, denned by the SPE e\ and €2, and 
the two-body part Hi, defined by the antisymmetrized, normalized tbme 
(ab\Hi\cd)j : T, of the Hamiltonian H. For later convenience, we also denote 
each matrix element by aij, bij, etc.. 





£1=0, 






e 2 = 10 




an = 


(111^111)0,1 = 


-3.0 


hi 


= <ll| J H>|ll> 1 ,o= 


-4.0 


ai2 = 


(ll|fl-j|12> ,i= 


-1.5 


bu 


= (111^/112)1,0=- 


-2.0 


013 = 


(ll|JTj|22) ,i= 


-1.7 


bi3 


= (ll|fl-j|22)i, = 


-1.9 


&22 = 


(12|ffj|12> ,i= 


-2.7 


b22 


= (12|^/|12)i, =- 


-2.9 


»23 = 


(12|#/|22) ,i= 


-1.6 


b23 


= (12|ffj|22)i >0 = 


-2.1 


«33 = 


(221^122)0,1= 


-2.5 


b33 


= (221^122)1,0= 


-2.9 


Cll = 


(121^112)0,0= 


0.0 


dn 


= (12|ffj|12)i,i= 


0.0 



Next we truncate the two-level space to a smaller space containing only 
the lower level, which means that we must construct the effective Hamilto- 
nian appropriate for this model space. 

As we mentioned in the Introduction, for the truncated space, the effec- 
tive interaction V e g for an ^4-body system can be written as 

v eS = v$ + V$ + ■■■ + . (10) 

(2) 

When only the two-body part V^ s is kept in shell-model calculations, the 
choice of the single-particle basis becomes important as we shall see below. 
In this section, we will use the original single-particle basis (|1) and |2), 
as used in Table 1), which is not a self-consistent basis, to calculate the 
effective two-body interaction V^ s ' . In the next section, we will introduce 
a self-consistent HF basis for A=A and re-calculate and compare the 
results obtained in these two sections. 

(2) 

To obtain the energy-independent effective two-body interaction V^ s , 
we start with the starting energy-dependent two-body G matrix and apply 
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Table 2: The A=2, 3, and 4 ground-state energies in the two-level model 
from exact matrix diagonalizations and from the one-level model space cal- 
culations with a two-body effective Hamiltonian obtained in the arbitrary 
and in the self-consistent HF bases. 



Ground-state energy 


Exact 


Arbitrary Basis 


Self-Consistent Basis 


E 2 (J=0, T=l) 


-3.390 


-3.390 


-3.390 


E 2 (J=1, T=0) 


-4.582 


-4.582 


-4.582 


£ 3 (J=l/2, T=l/2) 


-12.938 


-11.958 




E 4 (J=0, T=0) 


-27.703 


-23.916 


-27.645 



the Lee-Suzuki procedure (3|. Previously we labelled this as but now 
we relabel it as G JT to indicate the two-body conserved quantities. For the 
Hamiltonian and single-particle basis defined in Table 1, the G JT matrix 
element for the J=0, T=l channel for the chosen model space is given by 



G 01 (w 2 ) = an + ( ai2, a i3 ) 



^2 - (ei + e 2 + a 22 ) -a 2 3 \ / a 2 i 

-a 3 2 lo 2 - (2e 2 + 033) / V 031 



. / \ / u>2 -7.3 1.6 \ / -1.5 \ 
= -3.0 + -1.5, -1.7 ) 3 , „ • 1 

v / y 1.6 cj 2 - 17.5 y y -1.7 J 
Similarly, for the J=l, T=0 channel, we have 



G lu (w 2 ) = (-4.0)+ -2.0, -1.9 



cu 2 - 7.1 2.1 \ / -2.0 
2.1 uo 2 - 17.1 I \ -1.9 



(12) 



Clearly, these two-particle G matrix elements depend on the starting 
energy u>2- As in Ref.0, |[], we now use the Lee-Suzuki method || to obtain 
the energy-independent two-body effective interaction, which is equivalent 
to summing over the two-particle folded diagrams to all orders. We obtain 
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the following results: 




if lo 2 < 1.976 , 

if 1.976 < uj 2 < 12.595 , 

if uj 2 > 12.595 , 



(13) 



and 



(K (2), 10 



-4.582 if uj 2 < 1-288 , 

7.157 if 1.288 < to 2 < 12.391 , (14) 
17.625 if u 2 > 12.391 . 

It is easy to verify that the above numbers are precisely the eigenenergies 
of the two-particle system (A=2) whose Hamiltonian is defined by Table 1, 
as they must be ||, since it is just an application of Eqs.(^) and (0) to a 
simple A=2 case. 

The ground-state energies of the A=3 and A=4 systems, when we neglect 
effective many-body terms in Eq. (|l0|) , can be expressed in terms of these 
two matrix elements as 



3ei + 1.5 



o4 2) r + od 2) ) 10 



and 



E 4 = 4ei + 3.0 
The results are (noting that ei=0): 

£ 3 = 1.5(-3.390 - 4.582) 



-11.958 



and 



E A = 3.0(-3.390 - 4.582) = -23.916 



(15) 



(16) 



(17) 



(18) 



These differ from the exact results of -12.938 for E% and -27.703 for E4 by 
about 1.0 and 3.8, respectively. The discrepances reflect the importance of 
the neglected effective three- and (for A=4) four-body terms. 
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(2) 

4 Vg ff ; for a Self-Consistent Basis 

Addressing the case with the largest discrepancy, we now consider a self- 
consistent HF single-particle basis for the ground state of the A=4 system. 
We rewrite the original Hamiltonian of Table 1 as 

H = (H + U) + (Hj - U) = H'q + H'j (19) 

In the previous section, we have set U=0. In this section, we use a self- 
consistent one-body field U generated by all the A particles in the system. 
It can be obtained iteratively by using the following equations where n 
represents the iteration number: 

(2J + 1)(2T + 1) 



(a n \U n \(3 n ) = £ ^-p i /(l + ^,, 7 „_ 1 )(l + i_ 1 , 7 „_ 1 ) 

J,T,7=occ. \ Ja > 

{a. n -]_ln-\\Hl\^n-\ln-\) J,T , (20) 



(H + U n )\a n ) = e a Ja n ) , (21) 

where |aj), |/3j), |7j), etc. are the single-particle states in the i-th iteration 
and the summation (7) is over the occupied states. The tbme in the above 
expression are normalized and antisymmetrized. When converged, we obtain 
a self-consistent HF basis. 

It is obvious from the above equations that the resulting single-particle 
basis, defined by H' Q = (Hq + U), and the two-body residual interaction 
Hj = {Hi — U) are mass-dependent and should only be applied to the 
system for which they are calculated. 

It is easy to work out the self-consistent HF basis for A=4, which is a 
"closed-shell" system in our two-level model. The new single-particle states, 
denoted by |1') and |2'), are linear combinations of the old ones: 

1') = 0.92785|1) +0.37295|2), (22) 
|2') = -0.37295|1) +0.92785|2). (23) 
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The corresponding single-particle energies are: 

Hq\1') = -15.0323|1'), (24) 
Hq\2') = 5.3808|2'}. (25) 

The tbme of the residual interaction (Hj = Hj — U), evaluated using the 
new basis, are given in Table 3. 

Note that the tbme of Hj listed in Table 3 are ^-dependent not only 
because the one-body potential U is A-dependent but also due to the fact 
that when tbme are calculated for a one-body potential, a factor 1/(A — 1) 
has to be introduced: 

A 1 A 

i=l i<j 

With the self-consistent single-particle basis, there is no coupling be- 
tween the Op-Oh configuration and the lp-lh configuration. This is guaran- 
teed by the following equation: 

^2(2J + 1)(2T + 1)(1'1'|^|1 / 2 / )j iT = 0. (27) 
J,T 

Now in the HF basis, we can obtain an estimate of the ground-state energy 
of A=A even before calculating G JT by using the lowest-order in Hj estimate 
of y c ( ff 2) in Eq.flU): 

£HF = Ae v + 3 ((l'lW|l'l') 0il + (I'I'IaJII'iOlo) , (28) 

which gives £'| !F =-27.283 MeV. Note that this result is already closer to the 
exact energy of -27.703 MeV than the result obtained in the previous section 
for A=A. 

For the HF basis, we then calculate the effective two-body interaction 
for the truncated model space. The results are 

(V}g } ) 01 = 6.091 , (KffV = 4 -738 . (29) 
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Table 3: The single-particle energies of H' Q = (H + U) and the two-body ma- 
trix elements of Hi = (Hj — U) for the Hamiltonian defined in Table 1 with 
the self-consistent HF single-particle basis: [1'} = 0.92785|1) + 0.37295|2), 
and [2'} = -0.37295|1) +0.92785|2). 





£] 




-15.0323, 




£2 


= 5.3808 


(IT 


H'j 


IT 


o,i= 6.1415 


(IT 


\Hi 


1T)i,o= 


4.8074 


(IT 


H'j 


1'2' 


o,i= 0.0541 


(IT 


\Hi 


l'2')i,o= 


-0.0541 


(IT 


H'j 


2'2' 


o,i=-1.3402 


(IT 


\Hi 


2'2'}i, = 


-1.6120 


(1'2' 


Hi 


1'2' 


o,i= 4.5702 


(1'2' 


\Hi 


l'2')i,o= 


4.2265 


(1'2' 


Hi 


2'2' 


o,i= 1.2156 


(1'2' 


\Hi 


2'2')i, = 


0.8955 


(2'2' 


Hi 


2'2')o,i= °- 7399 


(2'2' 


\Hi 


2'2')i, = 


0.8176 


(1'2' 


Hi 


1'2' 


>o,o= 6.5505 


(IT 


\H'j 


l'2')i,i= 


6.5505 



So the ground-state energy for A=4, using an effective two-body interaction 
in a self-consistent basis, is 

E 4 = Ae v + 3.0 [{V$f l + (Vg*) 10 ] = -27.645. (30) 

This is extremely close to the exact result of -27.703 MeV. 

It should be pointed out that the good agreement with the exact result 
is mainly due to the fact that, with the HF basis, the matrix elements 
(lT|i^|T2') 0i i and (l'l'|J9 r £[l'2')i,o are very small (see Table 3). As we have 
mentioned previously, the sum of these two matrix elements vanishes in the 
HF basis [Eq.(^7])]. We further notice that the exact ground-state energy 
(i.e., the energy obtained in the full-space matrix diagonalization) of the 
A=A system, £4, does not depend the magnitude of each of the two matrix 
elements, which we denote by x: x= |(iT|F;|T2')o,i| = |(l'l'|.Hj|l'2 / )i,o|. 
However, it is clear that the effective 2-, 3-, and 4-body contributions to 
E4 do depend on x. In fact, as shown in Fig.l, the leading-order diagrams 
for the effective 2-, 3-, and 4-body forces involve x. Therefore, when we 
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increase x, the effective 2-, 3-, and 4-body interactions wiil ail change but 
these changes produce no net effect on E4 which remains the same. In 
particular, one can show that as x increases, the effective 2-body interaction 
becomes more attractive while the effective 3-body interaction becomes more 
repulsive. The smallness of x in the HF basis means that the effective many- 
body interactions are minimized and, consequently, the result using only the 
effective two-body interaction is in good agreement with the exact result. 

5 Conclusions 

With the Lee-Suzuki iterative method, it is now feasible to calculate the ex- 
act, energy-independent, effective two-body interaction for a no-inert-core 
model space from the starting-energy-dependent Bruckner G matrix G^(uj). 
We have demonstrated through a simple two-level model that when only the 
two-body effective interaction is used for shell-model calculations in a trun- 
cated model space, the choice of the single-particle basis is very important. 
An optimal choice of the single-particle basis should be chosen to minimize 
the neglected three- and more-body effective interactions. We have shown 
that a self-consistent HF basis serves this purpose very well for the two-level 
model. Encouraged by these results, it is now worth examining the more 
realistic situation when one must calculate G <y2 \uj2) in some chosen basis 
before any self-consistency calculation may be attempted. The present re- 
sults indicate that one should choose a realistic single-particle Hq for each 
A (assuming realistic is close to self-consistent) in order to minimize the 
neglected effective many-body forces. 
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Figure Caption 

Fig.l The leading-order diagrams for the effective 2-, 3-, and 4-body forces. 
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